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Abstract 

This work deals with a scalar nonlinear neutral delay differential equation 
issued from the study of wave propagation. A critical value of the coefficients 
is considered, where only few results are known. The difficulty follows from 
the fact that the spectrum of the linear operator is asymptotically closed 
to the imaginary axis. An analysis based on the energy method provides 
new results about the asymptotic stability of the constant and periodic solu- 
tions. A complete analysis of the stability diagram is given. Lastly, existence 
of periodic solutions is discussed, involving a Diophantine condition on the 
period. 
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1. Introduction 



In this work, we consider the following nonlinear neutral delay differential 
equation (NDDE) and we look for the solution y G Hl^j^{—1, +00), M), 



y'{t) + cy{t-l) + f{y{t)) + g{y{t - 1)) = s{t), t > 0, 
y{t) = yo{t) G H\{-1, 0),M), -1< t < 0, 

under the following assumptions: 

( c = ±l, 

feC\m,R), />o, /(o) = o, 

1^1 < 7I/I, with < 7 < 1, 
s{t + T) = s{t), T = — > 0. 

OJ 

In some cases, additional assumptions will be required: 

l^'l < 7l/'|, with < 7 < 1, 
lim fiy) = ±00. 



(1) 

(2a) 
(2b) 
(2c) 

(2d) 

(3a) 
(3b) 



Equation ([T]) is typically issued from hyperbolic partial differential equations 
with nonlinear boundary conditions 11, 0]- A closely related system of neu- 
tral equations has been derived, for instance, in the case of elastic wave 
propagation across two nonlinear cracks IJ] : y is then the dilatation of the 
crack, the shift 1 is the normalized travel time between the cracks, / and g 
denote the nonlinear contact law, and s is the T-periodic excitation. 

A natural issue when dealing with ([1]) is to prove existence, uniqueness 
and stability of periodic solutions. Articles and reference books [H, Ig], 0, 
19I, M, S, |16| usually address this question by considering a linear version of 
([T]) with |c| < 1. The critical value |c| = 1 raises technical difficulties: the 
spectrum of the linearized operator is asymptotically closed to the left of the 
imaginary axis, hence exponential stability is lost and small divisors may be 
encountered. 

An old contribution to the case |c| = 1 is done in j^. In the linear case, 
algebraic rate of convergence is proven. In the nonlinear case with c = — 1, 
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algebraic convergence is also proven in this reference, but assuming small 
data. In both cases, the main tools involved are asymptotic expansions of 
characteristic roots, Laplace transforms and function series. Another recent 
contribution for linear degenerate retarded differential systems is given in 
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The aim of our article is to push forwards the analysis of ([T]) in this 
critical case, especially to remove (when possible) the assumption of small 
solutions, and also to consider less regular initial data. A strategy based on 
energy analysis is followed. In the linear stability analysis shows that 

the energy method is optimal. When the stability condition is violated, the 
existence of an essential instability is proven ll|, [17 . 

The paper is organized as follows. Section [2] treats the full nonlinear 
NDDE. The stability of the zero solution in the homogeneous case is proven. 
Generalization to the non-homogeneous case with constant forcing is done. 
Non-constant periodic forcing is investigated in the case of small solutions. 
In section |3l new results are given in the linear case, concerning the stability 
diagram and the asymptotic rate of convergence. In sectionHJ the existence of 
periodic solutions is addressed. For particular values of the period, existence 
is proven whatever the amplitude of the source. More generally, existence is 
obtained under a Diophantine condition. This latter condition is satisfied if 
the period of the forcing is rational: T = p / q; if p is odd, we also exhibit a 
smoothing effect. Lastly, conclusion is drawn in section and some future 
lines of research are proposed. 



2. Nonlinear NDDE 

2.1. Stability in the homogeneous case 

Let us introduce basic functions and inequalities useful to state a main 
result: 

F{y)= r f{z)dz, G{y)= r g{z)dz, H{y) = 2 {F{y) - cG{y)) , 
Jo Jo 

(4) 

with /, g, 7 and c defined in ([I])-©- 

Lemma 1. For all y, F and H satisfy the inequality 

< 2 (1 - 7) Fiy) < H{y) < 2 (1 + 7) F{y). (5) 
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Proof. Let assume y > 0. From ([2]) and (jl]), it follows 

\G{y)\< r\g{z)\dz< [\\f{z)\dz= [\ f{z) dz = ^ F{y). 
Jo Jo Jo 

Similarly, if y < z < then \f{z) \ = —f{z) and 

\G{y)\< C\g{z)\dz< [\\f{z)\dz = - f{z)) dz = j F{y). 



Consequently, one gets < |G| < 7-F for all y. Using and c = ±1 
concludes the proof. □ 

Theorem 1 (Asymptotic stability of the zero solution). Let y be the solution 
of (OP without source term: s = 0. For all t > 0, one defines 

E{t) = (y\r) + /(y(r)))' dr > 0. (6) 

Then one obtains 

r+oc 

sup E{t) + sup F{y{t)) + / f{y{t)f dt < +oo, (7) 

t>0 t>0 Jo 

where F is defined in It follows the asymptotic stability of the origin: 

lim y{t) = 0. (8) 

t— S- + CXD 



The stability is proven by an energy method. Notice that the inequality 
d?]) allows convergence towards zero with an algebraic rate. Indeed there is 
not exponential stability of the equilibrium y = 0, as we will see for instance 
in the linear case in section [31 

Proof When s = 0, the NDDE writes 

y'{t) + f{y{t)) = -c (y'{t - 1) + cg{y{t - 1))) , (9) 
where c = ±1. Taking the square yields 

y\t) + = {y\t - 1) + cgiyit - 1)))' . (10) 
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Based on the right-hand side of f llUI) at t — 1 is rewritten 

{y' + cgY = {{y' + f)-{f-cg))\ 

= {y' + lf + {f-c9)'-2{y' + f){f-cg), 

= {y + f)' + if - cgf - 2y' if - eg) - 2 f if - eg) 



{y' + fy + {f-cg){f-eg-2f)-j^ H{y), 

[y^fy-[P-g^)-iH^y). 



The latter equation is injected into (fTOl) . which yields 



y\t) + f{y{t))) + {f{y{t - 1)) - 9\y{t - 1))) 
+|if(y(t-i)) = (y'(t-i) + /(y(t-i))'' 
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(11) 



(12) 



(13) 



Integrating (fT2|) over [0, t], we get for alH > 1 

£ + /(y(r)))' rfr + ' (f (y(r)) - /(y(r))) dr 

+H{y{t - 1)) - i/(y(-l)) = ' [y'{r) + /(i/(r)))' dr. 

Equation f|T3|) is simplified into the energy equality 

^ {y'{r) + /(y(r)))' + ' ifiyir)) - g\y{r))) dr + if(y(t - 1)) 
Z/o(r) + fiyoir)))' dr + if(yo(-l)) - r (f (yo(r)) - /(yo(r))) rfr. 



(14) 

From lemma [1] and equation (fT4|) , one deduces the energy inequality for all 
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T > 1 



1^'^ {y'ir) + /(y(r)))' dT + {l- 7^) 1^ ' fivir)) rfr + 2 (1 - 7) - 1)) 

< f {y'oir) + /(2/o(r)))' rfr + 2 (1 + 7) F(yo(-l)) - (1 - f ) j' f{yo{r)) dr 

(15) 

For alH > 1, it follows three inequalities: 

W))<^(^, ny{t))dt<^^, E{t)<C,. (16) 

The first inequality in f|T6|) implies y G L°°(0, +00). On any bounded set K, 
there exists d > such that |?/| < \ f{y) \ < |y| / don K. Using the second in- 
equality in flTBl) and the boundedness of y, one obtains y G i^^(0, +00). Lastly, 
the third inequality ensures that E is bounded; since f{y) G i^^(0, +00), one 
deduces that y is uniformly bounded in L^(T — 1, T), and hence ?/ is bounded 
in the Holder space C°'^/^(0, +00). It proves that y tends towards 0. □ 

Four remarks are raised by theorem [D 

• if c = and 7 = 0, then the NDDE is a simple ODE. Our computation 
provides for all T > 0: 2F{y{T)) + /(f (^'(t))^ + P{y{t))dt = 2F(y(0)) 
which implies the global stability of ?/ = 0. Moreover the stability of 
1/ = is exponential since /'(O) > 0; 

• in the usual case |c| < 1, then a stronger energy estimate can be ob- 
tained by our proof: y G H^iS^, +00)5 and the stability is exponential; 

• in the present case |c| = 1, then y E L"^ (1 L°° and y is uniformly 
bounded in H^{t — 1, r) for all r. As we will see in section [3] for the 
linear critical NDDE, the stability is not exponential; 

• if |c| = 1 and 7 = 1, then the stability is obtained, but the asymptotic 
stability can be lost. 
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2.2. Stability with constant source term 

In this subsection, we consider a non-homogeneous NDDE ([1]) with con- 
stant source term D: 

y'{t) + cy'{t-l) + f{y{t)) + g{y{t-l))=D, t>0, t^N, 

(17) 

yit) = yoit) e H\-l, 0), -1< t < 0, 

with assumptions (j2]). For the sake of clarity, we analyze successively g = 
and g 0. 

Corollary 1 (case g = 0). Let y be the solution of [T7\ ). Three cases occur: 

1. if D = f{d), then d is the unique globally attractive solution; 

2. if D > sup / (or D <'m.i f), then lim y{t) = +oo ( or — oo ); 

t->+oo 

3. if D = sup/ (or D = inf f), then there is no convergence towards a 
constant solution. 

Proof. We consider successively the three cases. 

Case 1: D = f{d). This case occurs if lim f{y) = ±oo. Injecting 

y— j-itoo 

z = y-d, f,^z) = f{d + z)- f{d) 
into f|T7|) yields the homogeneous NDDE 

z{t) + cz{t-l) + f',{z{t))=^- 



Since /d(0) = and > 0, the assumptions of theorem [T] are satisfied: 
lim z{t) = 0, and thus y tends asymptotically towards d. 

Case 2: D > sup f {oi D < inf /). To fix the minds, let us assume sup / < 
D < +00, and introduce 6 = D — sup / > 0. 

• If c = -1, then y'{t) - y [t - 1) = D - f{y) > 6, so that y' {n + r) > 
ynir) +n6, with r G [—1, 0]. As a consequence, lim y (t) = +oo, and 

hence lim y{t) = +oo. 

t— >' + 00 

• If c = +1, then y' (t) + y'{t — 1) > 5. Integration on [r — 1, r] yields 
?/(r) — y{r — 2) > 6, and once again lim y(t) = +oo. 

t— >-+oo 
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Case 3: D = sup / (or D = inf /). Let us assume y{t) d when t — )■ +00. 
Let yn{s) be y{n + s) for s G [0, 1]. We have ?/„(.) — )■ c? in L°°([0, 1]) and 
— in the sens of distribution. The NDDE can be rewritten as follows: 
y'^ + cy'^_i + fiyn) = D. Taking the weak limit, we get f{d) = D which 
is impossible since f < D. Consequently, y{.) cannot converge towards a 
constant. 

Using the same argument, we can state that y{.) cannot converge towards 
a periodic continuous solution. More precisely, let us assume that yn{s) — )■ 
d{s) where d{.) is continuous. Writing ?/„(0) = y{n) = ?/((n— 1) + 1) = 
yields d{0) = d{l) i.e. c? is a 1-periodic continuous function. Taking the weak 
limit in the NDDE yields a differential equation for d{.): (1 + c)d'(t) + 

fm) = D. 

• If c = —1 then f{d{t)) = D, which is impossible. 

• If c = 1 then d G C^([0, 1]) by the equation. Let Sq be a maximizer of 
d{.) on the compact set [0, 1], then d'{so) = and, by the differential 
equation, f{d{so)) = D which is again impossible. 

As a consequence, y{.) cannot converge towards a periodic continuous solu- 
tion. □ 



Now, let us examine the NDDE (|T7I) with all the terms. 

Corollary 2 (case g ^Q). Under the assumptions the unique d such as 
f{d) + g{d) = D is the unique globally attractive solution of [Tl\ ). 

Proof. Assumptions ([3]) easily imply that there exists a unique d satisfying 
f{d) + g{d) = D. Injecting 

z = y-d, fd{z) = f{d + z) - f{d), gd{z) = g{d + z) - g{d) 

into (fT7|) yields the homogeneous NDDE 

z'{t) + cz{t-i) + f'Mt)) + g'A^t - 1)) = 0. 

Inequality (I5al) implies that 7 > \gd\. Since /d(0) = and > 0, all the 
assumptions of theorem [T] are satisfied, and hence 2; — asymptotically. □ 

Without assumptions ([3]) in corollary [21 one may encounter more complex 
situations, with 2 or more solutions. Note that these assumptions are satisfied 



in the physically-relevant situation examined in [14 



8 



2.3. Stability with periodic source 

In this section, we consider the non-homogeneous NDDE ([1]) with the 
periodic source term s, but in the particular case (7 = 0. Existence of pe- 
riodic solutions to dl])-© is examined in section HJ Here, we assume that 
a periodic solution P{t) exists and the difference with another solution y{t) 
is bounded for all time . We focus on the convergence of y{.) towards -P(.)- 
The convergence is obtained under the condition ([Ti 



Proposition 1 (Conditional asymptotic stability of a periodic solution). Let 

y he the solution of (QP with f ^ and g = 0, and P be a periodic solution. 
Let I = conv +oo[, P(M)) be the convex hull of the two sets. If lis 

bounded and the condition 

^ei,deP{R)\f{w)- fid)-{w-d)f {d)\ ^ ' 

is satisfied, then lim {y{t) — P{t)) = 0. 

We try to extend the proof of Theorem [1] to more complex case, but the 
extension is very limited. It means that if P is a small periodic solution 
and y{.) is a small stable perturbation, then y{.) converges towards P. In 
other words, stability and condition ( fTSi) imply asymptotic stability. The 
difficulty follows from the occurrence of small divisors; see section [31 

Proof. Injecting z(t) = y(t) — P(t) in ([T]) yields 

z'{t)+cz'{t-l) + A{t,z{t)) = 0, (19) 

with 

A{t,z) = f{z+P{t))-f{P{t))= r f\u+P{t))du = z f f\sz+P{t))ds. 

Jo Jo 

(20) 

Let us introduce 



K{t, z)= f A(t, u)du. 
Jo 
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K is non-negative and convex with respect to 2;, and it satisfies 

I K{t, z) = z A{t, z) + P\t) [f\u + P{t)) - f{P{t))) du, 

= z A{t, z) + P'it) [f{z + Pit)) - f{P{t)) - z f\P{t))) , 

= z A{t, z) + P' (t) z'^ f"{zu + P{t)){l-u)du, 
= z A{t, z) + P it) z^ R{t, z). 

(21) 

Now, we mimic the proof of theorem [H Based on (fT^ and (PTl) . one deduces 

z\t) + A\t, z) + 2^ K{t, z)-2 P'{t) z\t) R{t, z)) = z\t - 1). (22) 

Since P is periodic and non constant, the term P' (t) z'^{t) R{t, z)) has no 
fixed sign and there is no hope to control (1221) . It motivates the condition 
fITS]) for small derivative. Based on fl2Ul) and fl^Tl) . this condition is equivalent 
to say that there exists < r < 1 such that 



2|P'(t)| 
^2\P'{t)\ 



f{z + P)- f{P) -zf\P) <T {f{z + P)- f{P)f , 

f" {uz + P{t)){l-u)du <t(^zJ^ f'{zu + P{t))dv^ 
^2P'{t) z\t)R{t, z)) < TA\t, z). 
Integrating (122|) over [0, r] and using the inequality ( l23l) yields 



(23) 



r-l 



z'\t) + (1 - r) / A\t, z)dt + 2 K{t, z{r)) 



< 



4'(t)c?t + 2K(0, Zo{0)). 



(24) 
(25) 



The first term in the left-hand side of (124|) implies that z is uniformly con- 
tinuous, whereas the third term implies that z is bounded. As in the proof 
of theorem [H the second term allows to conclude. □ 

The stability of the solution to ([1]) with a non-constant source term re- 
mains an open question. Such a L°° bound would imply that the interval / 
in proposition [T] is bounded. 
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3. Linear NDDE 



In this section, we consider the linear homogeneous NDDE issued from 

m 

y\t) + cy{t - 1) + ay{t) + hy{t - 1) = 0. (26) 

3.1. Stability diagram 

The stabihty of in f l26|l is analyzed in terms of a = /' and h = g . 
The cases |c| < 1 and |c| > 1 are well known in the literature, yielding 



to asymptotic stability and instability, respectively We focus therefore 
on the critical case |c| = 1 for linear scalar differential- difference equations 
discussed in [H, p. 191. In this book, it is said that the stability is too much 
intricate, so the model would never be used by engineers ! 

This analysis has two interests. First, it highlights the optimality of the 
energy analysis performed in theorem (T) as in the nonlinear case, where 
\g\ < f and / > 0, asymptotic stability is proven iff |6| < a. Second, the 
lost of exponential stability is explicitly exhibited. 

The stability analysis will involve the roots A of the characteristic equation 

h{\) = h{\; a,b) = X {e^ + c) + ae^ + b = (27) 
obtained by injecting a particular solution y{t) = e^^ into f l26|) . 

Proposition 2. The characteristic roots of satisfy the following prop- 
erties: 

1. the roots are located in a vertical strip Ainf < 3?(A) < Asup, where Ainf 
and Asup o'^e real constants depending on {a,b,c); 

2. the set of roots is countable infinite, yielding a sequence A„; 

3. for large values o/|A|, one gets the asymptotic value A„ ~ ifl{n), where 

f](n) = (2n + (1 + c)/2) tt; (28) 

4. the sign c)/3f?(A„) is constant on each quarter plans separated by a = ±b; 

5. the sign of 3?(A.„) changes when performing a symmetry relative to the 
line a + cb = 0. 
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Proof. We prove successively the properties. 

Properties 1 and 2. These are classical results 111, ll^, not repeated here. 
Property 3. If A is a root of f l27|) and |A| — +oo, then 



\ c \ + b , , 

e' = --— c, 29 

A + a 

and hence A = A„ ~ i Q{n) in ( 128|1 . 

Property 4. Let us consider a imaginary characteristic root X = iw. Injecting 
this root into (1271) leads to 

b + a cos w — w sinw = 0, (30a) 

a sinw + w (c + cos w) = 0. (30b) 

If sinw = 0, then w = nn and the following alternative occurs: if n = 
then a + b = 0, else a — cb = 0. If sinw 7^ 0, then one injects a sinw = 
—w (c + cos w) into ( l30i) . which yields 

6 sin w = —a sin w cos w + w sin^ w = cw cos w + w cos^ w + w sin^ w, 
= w {1 + c cos w) , 
= —ac sinw, 

(31) 

and hence a + cb = 0. By Rouche's theorem, the roots of depend 
continuously on a and 6. Consequently, the sign of 3fJ(A„) is constant on each 
quarter plans separated by the lines a = ±6. 
Property 5. From the characteristic equation fl27|) . one deduces 



h{X;a,b) = e-^ (cA (e"^ + c) + (6e~^ + a)) , 

= -ce^ {-X {e~^ + c) -cbe~^ -ca) , (32) 
= —ce^h{—X; —cb, —ca), 

which concludes the proof. □ 

For the sake of clarity, we will distinguish the cases c = +1 and c = — 1 
in (l26l) . The case c = +1 is examined in detail in the next theorem, whereas 
the modifications induced by c = — 1 are given further with a short proof. 
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Theorem 2 (Stability diagram for c = +1). We determine the stability of 
y = in IjiMB^ in terms of {a,b). Only one open quadrant is asymptotically 
stable, the other three being unstable: 

• a > \b\: asymptotic stability (in, 3?(A„,) < 0); 

• 6 < — |a|.- instability (in, 3ft(A„,) > Oj; 

• h> \a\: instability (1R(A„) > 0, except one root 3ft(A) < Oj; 

• a < —\h\: instability (^9fJ(A„) < 0, except one root 3?(A) > 0^. 
On the four edges, one has: 

• a = \b\: stability (1R(A.„) = 0, except one root A = —a < ifb>0); 

• —a = b> 0: stability (the spectrum belongs to iW, in, 9fJ(A„) = 0^; 

• a = b < 0: instability: 3ft(A„) = 0, except one positive root \ = a > 0. 
Finally a = 0, b = is a stable case. 

Proof. Based on property 4 of Proposition [21 we analyze successively the 
signs of 3f?(A) on each quarter plan. 

Quadrant 1: a > \b\. The proof of theorem [T] is adapted to the case of com- 
plex solutions y{t) = e^^, A G C. Instead of ( !T4|) . one obtains 

._l)|2 + (a2-62 
\yo{r)\'dr + {a-b)\yo{-l)\ 









f 


y\r) - 


■f ay(r) 









y{r)\'dT 





2 






dr + 





(33) 

This equality yields stability results when a — b > and — b"^ > 0, ie 
a > \b\. lfa> \b\, then asymptotic stability follows directly from fl33|) . as in 
the proof of theorem [TJ If a = \b\, the energy equality ( l33l) becomes 



y\r) + \b\y{r) dr + {\b\ - b) \y{t - l)f = Co- 



(34) 



Three cases are distinguished: 
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1. a = b = 0. The NDDE becomes y'{t) + y'{t - 1) = and h{s) = iff 
s(e'*' + l) = 0. The solutions of the characteristic equation (127|) are also 
known explicitly: A = and A = z (2 + 1) tt, /c G Z with multiplicity 
1, which proves the stability. Notice that the energy equality flM|l only 
means that all A„ G zM. 

2. a = 6 > 0. The NDDE becomes z'{t)+az{t) = with z{t) = y{t)+y{t- 
1). One therefore obtains z{t) = Ce~"* and y{t) = D~"'^+w(t), where w 
is a solution of the homogeneous difference equation w{t)+w{t — l) = 0. 
The spectrum is given by A = —a andA = i{2k + 1) n, k E'L. All the 
roots have multiplicity 1, which proves the stability. 

3. < a = — 6, energy yields 



which allows to conclude since y{t) is bounded. Furthermore, this 
equality also means that each characteristic root is simple and is an 
imaginary number. 

Quadrant 2: h < —\a\- Based on property 5 of proposition [21 the inner quad- 
rant and its two edges are successively analyzed: 

1. h < — |a|. From the case 1 of quadrant a > \b\, one deduces that 
5R(A) > when b < -\a\. 

2. b = a < 0. This line is obtained by symmetry of 6 = a > 0. From the 
case 3 of quadrant a > \b\, it follows that all the roots are imaginary, 
except one root A = —a > 0: this edge is therefore unstable. 

3. b = —a < 0. Symmetry of roots relative to a + b = yields 3fJ(A) = 0: 
this edge is therefore stable. 

Quadrant 3: b > \a\. When a = b = 0, the spectrum of ( 1271) is known: A = 
and A = Am = im = i{2k+l)7r. In the case A = 0, the implicit function 
theorem yields A (6) = —1/2 and hence A(6) = —6/2 + 0{b^). Consequently, 
one obtains 5R(A) < since 6 > 0. 

In the cases A = Am, the implicit function theorem yields A^(6) = —im 
and hence 




\m{b) =im b + Ab^ 



+ 0{b^). 



(35) 



To determine A, we inject ( 135|) into ( l27ll : 




(36) 
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Xmib) = tm--b+ {^^- — ) b^ + Oib^), (37) 



which leads to A = 1/(2 m^) — i/m? . Injecting A into fl35p yields 

- ^ + hr^ - - 

which therefore proves 3fJ(Am) > 0. 

This analysis of spectrum is valid only on a neighborhood of (a = 0, 6 = 0) 
with 6 > 0. The property 4 of proposition [2] allows to extend the result to 
the whole quadrant. 

Quadrant 4: a < — This case is obtained by performing a symmetry of 
roots on quadrant 3 relative to the line a + 6 = 0. It follows that there 
exists a single root with 3ft(A) > and an infinite sequence of roots such that 
3f?(A™) < 0. 

Edge b = —a > 0. Symmetry relative to a + 6 = implies that if A is a root of 
( 1271) . then —A is also a root. Moreover, is a simple root if 6 7^ 2; otherwise, 
h'{0, -6, 6) = 2 - 6 = and is of muhiplicity 2. 

Now, let us consider the points M(a, b) with — 6 < a < and O the 
origin. The segment OM is outside the half-line (—6, b) and belongs to the 
quadrant b > \a\. At M, h{0, a, b) = b — \a\ > and lim h{x, a, b) = — cxd: 

x— >■— 00 

there exists a negative real root Ao(a, b) at M. Since there is only one root 
with negative real part on quadrant 6 > |a| (see quadrant 3), Aq is unique. 
By symmetry, there exists also a unique positive real root in the quadrant 
a < —\b\. Consequently, the only root with negative real part on quadrant 
b > \a\ connects the only root with positive real part on quadrant a < — 16| 
at Ao(-fe, b) = 0. 

On the other hand, the other roots on quadrants b > \a\ and a < —\b\ 
satisfy 3?(An) > and 3?(A„) < 0, respectively. By Rouche's theorem, they 
depend continuously on (a, b). Consequently, 3ft(An) = on (—6, b). □ 

Four remarks are raised by theorem [2) 

1. The asymptotically stable region a > \b\ in the plane of parameters 
(a, b) is exactly the region given by theorem [H i.e. the energy method 
is optimal for the linear case to get the asymptotically stable region; 

2. a = \b\: the stability follows from the same arguments as in the case 
a > \b\: energy method: theorem (H and a study of the characteristic 



roots as in 



14l ]. Moreover, crossing this edge, all real parts become 
positive in b < —\a\. This kind of dramatic transition, where all real 
parts are negative on a > \b\ becomes positive on 6 < —\a\, is an 
essential instability (ill . 0, flTj : 
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3. —a = b > 0: this is a stable case if and only if all roots are simple. But 
it is not always the case. For instance, the root s = is simple on this 
edge except (a = —2, b = 2). In this last case, s = is of multiplicity 
2 and the unbounded solution y{t) = t occurs; 

4. the central case a = 6 = 0, at the boundary of all previous subdomains, 
is stable. 

We treat briefly in a similar way the case c = — 1. Stability of y = in 
is studied in terms of (a, b). 

Theorem 3 (Stability diagram for c = —1). 

Only one open quadrant is asymptotically stable, the other are unstable: 

• a> \b\: asymptotic stability (in, 3ft(A„,) < Q); 

• b> \a\: instability (Vn, 9fJ(A„) > 0); 

• b < —\a\: instability (^3fJ(A„) > 0, except one root 5R(A) < 0); 

• —a > \b\: instability (^{Xn) < 0, except one root 5R(A) > 0). 
On the four edges, one has: 

• a = \b\: stability (1R(A„) = 0, except one root A = —a < ifb<0); 

• a = b> 0: stability (1R(A„) = 0); 

• a = —b < 0: instability (1R(A„) = 0, except one positive root A = —a). 
At the origin, one has 

• a = b = 0: weak instability (^{\n) = 0, except with multiplicity 2). 



Proof. The proofs of the various cases are shortly sketched: 
1. We begin by the energy equality for c = —1: 









2 




/' 




f ay(r) 


dr -(- 


-ia + b) 1 













t-l 



\y{r)\'dT 



y'o{r) + ayo{r) dr + / \yo{T)\^ dr + {a + b) \yo{-l)\' 



This energy equality yields the asymptotic stability when a > \b\. 



(38) 
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2. When a = -b, the NDDE becomes z'{t) + az{t) = with z{t) = 
y{t) — y(t — 1). The spectrum is therefore —a and on the imaginary 
axis. When a = 0, is a root with multiphcity 2; when a ^ 0, the 
roots are simple. 

3. The stabihty on a = 6 > is given by the energy equahty. 

4. The symmetry induced by property 5 of proposition [2] concludes the 
case 6 > |a| and b = —a > 0. 

5. Passing trough a = 1,6 = —1, we obtain the sign of the real part of 
roots in 6 < —\a\ for instance, and by symmetry in —a > \b\. 

6. On a = b < 0, the real part of roots are null, as in the proof of theorem 

13 

□ 

3.2. Small divisors 

Now, we consider the linearized non- homogeneous NDDE issued from ([T]) 

Cy = y\t) + cy'it - 1) + ay{t) + by{t - 1) = (39) 

where s is a T-periodic function or an almost periodic function jsj, and 
c = ±1. We assume a > \b\: as stated in theorems [2] and |3l this is the only 
asymptotically stable case. Based on ( 127|) . we introduce "H = e~'^ h{X) which 
satisfies 

£(exp(At)) = n{X) exp(At), 

niX) = A (1 + cexp(-A)) +a + b exp(-A)), 

and iTii^iu)] = \h{iu)\. Lastly, we denote by S and Y the mean values of s 
and y over one period: for instance, Y{t) = Jt-Tl/i''') 

If S* = 0, then elementary calculations yield CY = 0. In this case, 
theorem [1] ensures that lim Y{t) = 0. The mean value of the periodic 

solution - if it exists - is therefore asymptotically stable. 

If the source is monochromatic s{t) = e*"^*, then y{t) = yp{t) + z{t), 
where the unique periodic solution is Vpit) = e*'^* / T-L^iu), and z satisfies the 
homogeneous equation (126|) : Cz = 0. Theorem [1] ensures that lim z(t) = 0. 

In other words, a monochromatic source provides a unique periodic solution 
which is globally asymptotically stable. 

Now, we consider a source with an infinite spectrum: 

s(t) = 5^Sfce^-^*. (41) 
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A formal solution of fl39D is 



fcez 



The denominator in (142|) never vanishes on the imaginary axis {H{iuj) 7^ 0), 
but 

liminf lH^iu)] = 0, limsup iTii^iu)] = +00, 

and hence various cases can occur. We illustrate some situations on particular 
values of the period of the forcing T: 



if T = 1, then Uk) = 2/c7r) = 4i kir + a + b, and hence 
|'H(za;fc)| — +00 as /c — )■ +00. Consequently, a smoothing effect 
occurs: s E L"^ leads to i/p E H^; 

HT = 2, then n{iuJk) = n{ik7r) = ikir (l + (-1)'=) +a+ (-l)H 
and hence T-L^iUk) = a — b if k is odd: no smoothing effect. 



— if T = 1, then T-iiiuk) = a + b: no smoothing effect; 

- if T = 2, then?/(ia;fc) =ikiT (l - (-1)'=) +a+ (-l)H and hence 
l-LiiuJk) = a + 6 if is even: no smoothing effect. 

A worst behavior is expected if lim T-Lliuk) = 0, where small divisors 

may occur. The goal of the next theorem is to determine the sequence of 
frequencies Uk leading to this situation. 

Theorem 4. We assume \c\ = 1 and a > \b\ in (2l\ ). If lim h{iuk) = 0, 
then there exists a function </> : Z — )■ Z and a real sequence Sk such that 

lim = +00, lim Sk = 0, 

where fl{k) = {2k + (1 + c)/2)7r is defined in ^2^) . and 

d = a-cb>0. (44) 
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a2-62 / 1 \^ 

h{lUk) = -Sk + i^rTrT-Tj^ +0{\Sk\ + T^rT-Tj--T] . (45) 



Moreover, the small divisor has the same modulus than 



Proof. Two steps are involved: 

Step 1: h{iujk) 7^ 0. Under the hypothesis a > \b\, the property 1 of propo- 
sition |2] can be made more precise. The characteristic roots f l27|) satisfy 

(A + a)e^ = -(cA + 6). (46) 

The equahty A + a = imphes b = ca which is impossible when \b\ < a. As 
a consequence, one obtains 

e^ = -^^ Kl= <1, (47) 
A + a ' ' 

which proves that 3ft(A) < 0. The characteristic function h therefore never 
vanishes on the imaginary axis. 

Step 2: asymptotic expansion of Uk- From (!27|) . it follows 

l + ce-- + of^V (48) 



By hypothesis, h{iuk) — when \uk\ — +oo. It results that 1 + ce — 0, 
and hence there exists Zk = Q{4>{k)) and — > so that 

uJk = Zk + rk, (49) 



with flk defined in (j28|) . Injecting P9|) into fl271) and using e^*^'' = — c, we 
obtain 

hitcok) = t{z, + rk) {l-e~'"^)+a-cbe-'"', 

= i{zk + rk) (^irk + '^ + Oirl)^ +a-cb {l-irk + 0{rl)) , 
= -ZkTk + icbrk + a - cb + O [z^rl] + O (r^) . 

(50) 

Since h{iuk) — ?■ 0, there exists Sk such that 

rk = \ , (51) 

Zk Zk 
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with d defined in (jUj). Putting (15T]) into the second hne of (13U]) gives 

M...) = + A + cm) + o (I) + o (i^) + o(4)^ 

Notice that (f/2 + chd = d{d/2 + cb) = d {{a- cb)/2 + cb) = d{a + cb)/2 = 
(a — c 6) (a + c 6)/2 = (a^ — 6^)/2 > 0, which concludes the proof. □ 

Theorem H] states that the small divisors problem occurs when the spec- 
trum of the source s approaches the asymptotic values of the characteristic 
roots in the way fH3|) . If s is T-periodic, then the solution of f l39|) looses 
accordingly at most one derivative. 

4. Existence of periodic solutions 

We consider again the full nonlinear non autonomous NDDE ([T]): indeed, 
some results for the nonlinear case are based on the linear analysis performed 
along section [31 We denote by the Sobolev space of T-periodic functions 
with derivatives in L^. 

Two cases are distinguished, depending on the period T of the source. If 
1 / T G N, then general results can be proved whatever the amplitude of the 
solution. For a much larger set of periods and for small sources, results can 
be obtained under a condition on the pulsation u = This condition is 
equivalent to a Diophantine condition on T. We finish the section by some 
remarks and open problems about this Diophantine condition. 

4.I. Periods T = 1 / n 

Theorem 5. Properties ^ and Iji^ are assumed. If nT = 1, n & N, and 
the source s{.) belongs to H^, then there exists a unique T-periodic solution 
y G to (Cy, which satisfies the non autonomous ODE 

(1 + c) y\t) + fiy{t)) + giy{t)) = sit). (52) 
Furthermore, if c = 1, we can take the source G Lj. 

In other words, if the delay 1 is a period of the source, there exists only 
one periodic solution for the (NDDE) with the same period as the source. 
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Proof. Considering a 1-periodic solution y to ([T]) gives fl52p . If c = +1, then 



y'(t) + %(t)) = ^, (53) 

with h = (/ + (?)/ 2. Under the assumptions ([3]), h satisfies h{0) = and 
h (0) > 0. The existence, uniqueness and asymptotic stabihty of a periodic 
solution to (|53|) follows accordingly; see for instance the proof of theorem 5-1 
in [13] • If c = — 1, then ([3]) ensures that s belongs to the range of h, and 
hence y{t) = h-^{s{t)). □ 

Theorem [5] does not ensure the asymptotic stability of the solution to 
the NDDE ([T]). For this purpose, an hypothesis of stability of small solution 
must be added, as in proposition [TJ 

4-2. General periods 

Unlike the special cases investigated in section 14. 1^ the results obtained 
in this section involve a much larger set of periods T but for small sources. 
The next theorem is also valid for almost periodic sources: it suffices to 
replace the periodic spectrum {ku, A; G Z} by the almost periodic spectrum 
{uk, k E Z} in (l54l) . ( 155|) below. Quasi-periodicity can occur, for instance, 
when the source is made of multiple incommensurable periods. Now, we state 
the main result of existence. 

Theorem 6. Let us consider the nonlinear non autonomous NDDE (J\) with 
assumptions ^-^3^). Based on the integer Qk ^ME), the subsequence (j) is 
chosen such that 

\ku-n{(j){k))\=min\kuj-n{n)\. (54) 

n 

The source is T -periodic, small and belongs to . If the condition 

lim inf |rf - fi(0(A;)) x{ku- fi(0(A;))) | > (55) 

k 

is satisfied, then there exists a unique periodic solution y G H^. 

Proof. Theorem m and condition fl55l) imply that infk \ h{iuJk)\ > 0. The 
operator s{t) — )■ y{t) in the linearized NDDE fl26l) is thus continuous. The 
implicit function theorem can therefore be applied, ensuring the existence of 



a unique solution in a neighborhood of the origin as in [ij]. □ 
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Note that (p^k) ~ |; — ?■ ±00. The condition (l55l) means that there are no 
small divisors in the linearized NDDE ( l39ll . Now, the key issue is to know 
when fl56|) or fl55l) are fulfilled. An important case is obtained for rational 
values of the period, as stated in the next theorem. 

Theorem 7 (Rational period). The same notations and assumptions are 
used as in theorem\^ If T & Q, then the condition 153]) is fulfilled and the 
conclusion of theorem holds. 

Furthermore, if c = 1, let us write T = p / q with gcd{p, q) = 1. If p is 
odd, then a small source in yields a unique periodic solution y G H^. 

In other words, no small divisors occur when the period of the source is 
rational. Moreover, when c = +1, periods in the form T = ^l^±l (m, q G 
Z, q ^ 0) provide a smoothing effect. 

Proof. Let Vk be min\kuj — Q{n)\ = \kuj — Q{(f){k))\. The condition (155|1 is 

nGZ 

fulfilled if there exists 6 > such that Zk = \d — Q{(j){k)) Vk\ > 6 > for \k\ 
large enough. Before examining the existence of 5, we recall that d > and 
n{(j){k)) ~ ku. 

The distance to a real subset A C M is denoted: distfx. A) = inf |x — a|. 



From the definitions of tu ( I2dl) and ( l28l) . it follows: 

k 1 + c 



Vfc = 2 vr min 



T 4 



n 



= 2 7rdist I — — ^— ^, Z , , 
\T 4 ' ' ' 

= 2 7rdist ( A;- - Z 1 . 

\ p 4 ; 

Since dist(xZ) < ^ for all x, it follows < t'jfc < vr. Moreover, periodicity 
of (wfc) implies that the sequence (f^) takes at most p values. Let < v = 
inf Vk = min v^. If v = 0, then we take 6 = d. If t; > 0, then any positive 5 

works since — t- +00. Consequently, the condition is always satisfied 
when T is rational, ensuring the existence of a unique periodic solution. 

Let us now examine deeply the cases where v_ = 0. This equality occurs 
in two cases: 

• If c = —1, then v = 0. It follows from Vkxp = 2% dist (k q, Z) = 0. 
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• If c = +1, then = if and only if p G 2Z. This necessary and 
sufficient condition follows from Vk = 2 7rdist ^A;| — |, = for 

some k ^ 0. It amounts to say that k ^ — ^ = m for some m G Z, i.e. 
q / p = (2 m + 1) / {2 k), and hence p has to be even. 

It proves that ?; = iff c = +1 and p is odd. In this case, the characteristic 
function n{X) in gO]) satisfies n{X) = \ {l + c e"^) + (a + 6 e~^) ~ A (1 + e~^) 
when \k\ — ?■ +oo. Since the linearized operator is a Fourier multiplier by the 
characteristic function, \'H{Xk)\ ~ \k\ ensures that the linearized equation is 
smoothing from into H^. □ 

4-3. About the Diophantine condition 

Up to now, we have proven that rational periods satisfy the condition 
( 156|) . ensuring the validity of theorem [61 Are there other sets of periods in 
a similar case ? We have not resolved this question. Here, we only propose 
few results on this topic. More precisely, we study the complementary of this 
set, when the small divisors problem occurs. To estimate this complementary 
set, the following definition is introduced. 

Definition 1 (Diophantine condition V). Let D = d / (4-7r^) > 0, with d = 
f (0) — eg (0) > 0. The real number 9 G satisfies the condition D if and 
only if there exists an infinite number of integers (n, k) such that 

The next lemma states that the condition (IS^ is exactly the converse of 
the condition (ISSj) . 

Lemma 2. The period T satisfies ([5^ if and only if 6 = ^ does not satisfy 

(Eg;. 

Proof. To prove the equivalence, let us rewrite the condition in theorem |4] 

to have a small divisors problem. In this case, u is such that there exists 
: Z I— >■ Z satisfying 

\<p{k)\ oc, 

kuj-VL{<p{k)) 0, 

Vt{(t){k)){kuj -Vt{(t){k)) d, 
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when \k\ — )■ oo. This is equivalent to say that there exists an infinite number 
of (/c, n) G Z X Z such that 



T 



1 



n 1+c DT 



k Ak F 




) 



(57) 



with D = consequently, 1 / T does not satisfy the Diophantine condition 



The size of the set of numbers satisfying the Diophantine condition (1561) 
is unknown. To advance on this question, we consider C G M, and we define 
the set E{C) as follows: 9 G E{C) if there exists an infinite number of 
(A;, n) G Z X Z such that 



E{C) is close to the set of reals satisfying Indeed, f +^ = ^^^^^ G Q 
and DT is replaced by C in f l57|l . Moreover, E{C) is bigger than the set of 
numbers satisfying V, since the condition fISBl) is more restrictive than fISS]) . 

Proposition 3. Depending on the value of C, we get the following alterna- 
tive: 

• if C = 0, then Q C E{C). Moreover, almost all real numbers belong to 
E{C), but E{C) ^ R; 

• ifC ^0, thenQnE{C) = (lS. 

Proof. We examine successively the two cases: 

• Case C = 0. Clearly all rational numbers are in -E'(O). 

We recall a Khintchine's theorem on Diophantine approximation (l5l.[5|. 
If is a non- increasing function, : N* H-jO, +oo[ and if ^„>o0('^) = 
+00 then, for almost all real numbers, there are infinitely many rational 
p / q such that \x — p / q\ < (j){q) / q. 

It suffices to take 0(n — 1) = 1 / (nln(n)) to conclude that almost all 
real numbers are in -E'(O). Indeed, we have 



dSSD. 



□ 




(58) 



,2 ij^q 



1 



= 




) 
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Notice that -E'(O) 7^ M. For instance, the golden ratio (1 + / 2 does 
not belong to E{0), see [15 . 



• Case C ^ 0: let 9 = p / q he in E{C). From (ISHD, it follows that 
kp — nq = ^ + o{k~^), and hence kp — nq is a. convergent sequence 
of integers towards 0. Such a sequence of integers is equal to for 
k large enough. But this is impossible since C 7^ 0. Consequently, 
Qn E{C) = 0. 

□ 

5. Conclusion 

The main results of this article are summed up as follows. 

• theorem [H corollary [1] and corollary |2) asymptotic stability of constant 
solutions in the nonlinear case is proven for a constant forcing. The 
results are valid whatever the amplitude of the source; 

• proposition [D asymptotic stability of a periodic source under a periodic 
non-constant forcing and in the nonlinear case is obtained under two 
assumptions (stability and small amplitude of the source); 

• theorems [2] and [3l stability diagram is analyzed whatever the coeffi- 
cients of the linearized NDDE; 

• theorem (5) existence of a periodic solution to the nonlinear NDDE is 
proven whatever the amplitude of the source but for particular values 
of the period. Theorem [6] considers a much larger set of periods, but 
under the assumption of small source and a Diophantine condition. 
Theorem [7| states that rational periods satisfy this condition. 

Deeper analysis is required about the set of real numbers not satisfying fl56|) . 
Proposition [3] is a first step in this direction, but the measure of the set 
involved remains an open question. 
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